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Method of evaluating a physical quantity represen-bative 
of an interaction between a wave and an obstacle 

The invention relates to the modelling of the 
5 interactions between an incident wave and an obstacle 
of this wave, in particular in the area of 
nondestructive testing. 

The modelling of the interactions between a wave and an 
10 obstacle receiving this wave, such as a target placed 
in the responsive zone of a sensor, finds an 
advantageous application in nondestructive testing. 

A method of modelling called "finite elements" is known 
15 consisting in applying a tiling of the three- 
dimensional space surrounding the obstacle and in 
evaluating the aforesaid interactions for all the tiles 
of the space. 

20 Methods of computation by "finite elements" afford a 
solution to a problem posed in the form of partial 
differential equations. They are based on a 
representation of the space under study by an 
assemblage of finite elements, inside which are defined 

25 approximation functions determined in terms of nodal 
values of the physical quantity sought. The continuous 
physical problem therefore becomes a discrete finite 
element problem where the nodal quantities are the new 
unknowns. Such methods therefore seek to approximate 

30 the global solution, rather than the starting equations 
in the partial spatial derivatives. 

The discretization of the space taken into account 
ensures that the latter is entirely covered by finite 
35 elements (lines, surfaces or volumes) , this operation 
is called "meshing" in two dimensional space (2D) or 
"tiling" in three-dimensional space (3D) . The elements 
involved are either rectangular or triangular in 2D, or 
parallelepipedal or tetrahedral in 3D. They may be of 
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different sizes, distributed uniformly or otherwise 
over the surface- 



In general, the physical quantity sought, such as an 
5 electrostatic potential or a pressure value, is known 
on the boundary of the domain. This boundary may be 
fictitious. Boundary conditions are imposed there. The 
potential is therefore unknown inside the same domain. 
A node is then defined as being a vertex of an element. 
10 The unknowns of the problem are therefore the values of 
the potential at each node of the domain as a whole. 

By way of illustration. Figure 6 of the prior art 
represents an exemplary surface, consisting of two 
15 materials Ml and M2, of different electromagnetic 
properties, and meshed by triangular elements each 
comprising three nodes Ai, Bi and Ci. The domain as a 
whole is delimited by a boundary F. 

20 Once the mesh has been defined, several approaches 
exist for transforming the physical formulation of the 
problem into a discrete modelling by finite elements. 
If the problem is formulated through differential 
equations and consists in minimizing a functional, then 

25 a variational procedure is generally applied. This 
transformation leads to a matrix formulation which when 
solved gives the nodal solutions, the solutions at the 
non-nodal points being obtained by linear 
interpolation . 

30 

Nevertheless, such computations, in three dimensions, 
require considerable computing resources and generate 
very long computation times, despite the enhancement in 
the performance of software allowing the implementation 
35 of these computations. 



Admittedly, 2D problems, often simplified by symmetry 
conditions that are advantageous for modelling only 
part of the geometry, are solved rapidly. However, this 



- 3 - 

is not so for 3D problems, which are the most frequent. 
Figure 6 shows how the fineness of the mesh, that is to 
say the ratio of the size of an element to that of the 
smallest detail of the domain, amplifies the number of 
5 nodes. 



Consequently, the number of equations and of unknowns 
increases proportionally, and, hence, the computation 
time required for solving the problem. It is important 
10 to point out that the generation of the mesh, namely 
the discretization of the workspace, and the generation 
of the list of nodes consumes greater computation time 
than that required for solving the problem. 

15 The present invention aims to improve the situation. 



Accordingly it proposes a method of evaluating a 
physical quantity associated with an interaction 
between a wave and an obstacle, in a region of three- 
20 dimensional space, in which: 

a) a plurality of surface samples, of which a part at 
least represents the surface of an obstacle 
receiving a main wave and emitting, in response, a 
secondary wave, is determined by meshing, and at 

25 least one source emitting an elementary wave 

representing a contribution to the said secondary 
wave is allocated to each surface sample, 

b) a matrix system is formed, comprising: 

- an invertible interaction matrix, applied to a 
30 given region of space and comprising a number of 

columns corresponding to a total number of 
sources , 

- a first column matrix, each coefficient of which 
is associated with a source and characterizes the 

35 elementary wave that it emits, 

- and a second column matrix, which is obtained by 
multiplication of the first column matrix by the 
interaction matrix and the coefficients of which 
are values of a physical quantity representative 
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of the wave emitted by the set of sources in the 
said given region, 

c) to estimate the coefficients of the first column 
matrix, chosen values of physical quantity are 

5 assigned to predetermined points, each associated 

with a surface sample, the said chosen values being 
placed in the second column matrix, and this second 
column matrix is multiplied by the inverse of the 
interaction matrix applied to the said predetermined 
10 points, 

d) to evaluate the said physical quantity representing 
the wave emitted by the set of sources in a given 
region of three-dimensional space, the interaction 
matrix is applied to the said given region and this 

15 interaction matrix is multiplied by the first column 

matrix comprising the coefficients estimated in step 
c) • 

Thus, according to one of the advantages afforded by 
20 the present invention, the meshing step a) relates only 
to one or more surfaces, whereas the method of 
modelling of the "finite element" type requires a 
tiling of the whole space neighbouring the obstacle, 
thereby making it possibly to reduce, in the 
25 implementation of the method according to the 
invention, the memory resources and the computation 
times required. 

The method according to the invention applies equally 
30 well to a main wave emitted by a far source as to a 
main wave emitted in the near field. 

Advantageously, to evaluate a physical quantity 
representative of an interaction between an element 
35 radiating a main wave and an obstacle receiving this 
main wave, 

- in step a) , a plurality of surface samples together 
representing an active surface of the element 
radiating the main wave is furthermore determined, by 
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meshing, and at least one source emitting an 
elementary wave representing a contribution to the 
said main wave is allocated to each sample of the 
active surface, 

5 - steps b) , c) and d) are applied to the samples of the 
active surface, and 
- the said physical quantity representing the 
interaction between the radiating element and the 
obstacle in a given region of three-dimensional space 
10 is evaluated by taking account of the contribution, 

in the said given region, of the main wave emitted by 
the set of sources of the active surface and the 
contribution of the secondary wave emitted by the set 
of sources of the surface of the obstacles. 

15 

The terms 
either an 
generator, 
sensor of 

20 

In a first embodiment, the physical quantity to be 
evaluated is a scalar quantity and, in step a) , a 
single source is allocated to each surface sample. 

25 In a second embodiment, the physical quantity to be 
evaluated is a vector quantity expressed by its three 
coordinates in three-dimensional space, and three 
sources are allocated, in step a) , to each surface 
sample . 

30 

In an advantageous embodiment, to estimate, in step d) , 
the contribution of the secondary wave in the given 
region of space, the values of physical quantity chosen 
in step c) are dependent on a predetermined coefficient 
35 of reflection and/or of transmission of the main wave 
by each surface sample of the obstacle. 

Thus, it will be understood that the secondary wave may 
either correspond to a reflection of the main wave, or 



"radiating element" are understood to mean 

emitter of the main wave, such as a wave 

or a receiver of the main wave, such as a 
this wave. 
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to a transmission of the main wave, or else to a 
diffraction of the main wave. In this advantageous 
embodiment, step c) corresponds finally to a 
determination of the boundary conditions at the surface 
5 of the obstacle, in the guise of interface between two 
distinct media in particular in a heterostructure . 

Furthermore, for nondestructive testing of a target 
forming an obstacle of a main wave, a chosen 

10 coefficient of reflection or of transmission is 
allocated to all the predetermined points of the 
surface of the target, and a simulation obtained by the 
implementation of the method within the meaning of the 
invention is compared with an experimental measurement, 

15 Thus, the points of the surface of the target which, in 
the experimental measurement, do not satisfy the 
simulation correspond to inhomogeneities or to 
impurities on the surface of the target. 

20 In another approach, the global properties of the 
obstacle are known, in particular in transmission 
and/or in reflection. By the implementation of the 
method of the invention, the position in space of a 
sensor or even the shape of this sensor is then 

25 optimized for application to nondestructive testing, 
this sensor being intended to analyse a target forming 
an obstacle of the main wave. 

Accordingly, in an advantageous embodiment, a plurality 
30 of values of the physical quantity estimated in step d) 
within the meaning of the method of the invention, 
which are obtained for a plurality of regions of space, 
are compared so as to select a candidate region for the 
placement of a radiating element intended to interact 
35 with the obstacle. 

As indicated hereinabove, the terms "radiating element" 
are understood to mean either a sensor or a generator 
of the wave. It will thus be understood that the 
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optimization of the position of the radiating element 
can be applied also to the optimization of the 
placement or of the shape of a wave generator. For 
example, the present invention also finds an 
5 advantageous application to the placement of 
loudspeakers in a closed volume, delimited by 
obstacles, such as for example the cabin of a motor 
vehicle . 

10 Other characteristics and advantages of the invention 
will become apparent on examining the detailed 
description hereinbelow, and the appended drawings in 
which: 

15 - Figure lA diagrammatically represents the respective 
surfaces of a radiating element ER emitting a wave 
and of an obstacle OBS receiving this wave, said 
surfaces being meshed with a view to evaluating a 
scalar quantity representative of the wave at a point 

20 M of three-dimensional space; 

- Figure IB represents in detail a surface sample dSi 
corresponding to a mesh cell of Figure lA, as well as 
a source Si associated with the surface sample dSi; 

25 

- Figure 2A diagrammatically represents the respective 
surfaces of a radiating element ER emitting a wave 
and of an obstacle OBS receiving this wave, said 
surfaces being meshed with a view to evaluating a 

30 vector quantity representative of the wave at a point 

M of three-dimensional space; 

- Figure 2B represents in detail a surface sample dSi 
corresponding to a mesh cell of Figure 2A, as well as 

35 three associated sources SAi, SBi and SCi; 

- Figure 2C represents, viewed face on, a meshed 
surface, each surface sample of which comprises three 



sources SAi, SBi and SCi, for the estimation of a 
vector quantity; 

Figure 3A represents, by way of illustration, the 
plates of a capacitor, of respective electric 
potentials VI and V2, for the estimation of an 
electric potential at the point M of three- 
dimensional space, a single source Si being 
associated with each surface sample dSi of Figure 3A; 

Figure 3B represents, by way of illustration, the 
plates of a capacitor, of respective electric 
potentials VI and V2 , for the estimation of an 
electric field E(M), at the point M of three- 
dimensional space, three sources SAi, SBi and SCi 
being associated with each surface sample dSi of 
Figure 3B; 

Figure 4A represents, like Figures lA and 2A, an 
interaction between a radiating element ER and an 
obstacle OBS, so as to evaluate a physical quantity 
(scalar or vector) at a point M in a portion of space 
delimited by the surface of the radiating element and 
the surface of the obstacle, this point of space M 
receiving both the wave emitted by the radiating 
element and the wave reflected by the obstacle; 

Figure 4B, complementary to Figure 4A, represents a 
transmission by the obstacle OBS of the wave emitted 
by the radiating element ER, at a point M of a half- 
space delimited by the plane formed by the surface of 
the obstacle OBS; 

Figure 5A diagrammatically represents an obstacle 
OBS, of finite dimensions, with sources associated 
with the surface samples and arranged so as to 
estimate a quantity representative of a reflection of 
the wave off the obstacle; 
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- Figure 5B, as a supplement to Figure 5A, 
diagraminatically represents an obstacle OBS, of 
finite dimensions, with the sources associated with 
the surface samples and placed so as to estimate a 

5 quantity representative of the transmission of the 
wave by the obstacle; 

- Figure 5C represents a simulation of an ultrasound 
wave emitted by a radiating element ER and 

10 propagating towards an obstacle OBS; 

- Figure 6 represents a mesh of three-dimensional 
media, for the application of a method of computation 
by "finite elements", within the meaning of the state 

15 of the art; 



- Figure 7A represents in detail a surface element and 
an observation point M whose relative positions are 
labelled by an angle 9; and 

20 

- Figure 7B diagrammatically represents a surface to be 
meshed of complex shape, in particular with an 
observation point M lying in a shadow zone with 
respect to certain sources of the surface. 

25 

Reference is first made to Figure lA, in which the 
surface of an obstacle OBS, receiving a wave, is meshed 
according to a plurality of surface samples dSi to dS4, 
in accordance with the aforesaid step a) • 

30 

Referring to Figure IB, with each surface sample dSi is 
associated a hemisphere HEMi, tangential to the surface 
sample dSi at a point of contact Pi. Preferably, this 
point of contact Pi corresponds to the apex of the 
35 hemisphere HEMi. For the estimation of a scalar 
physical quantity at the point M (such as an 
electrostatic potential, an acoustic pressure or the 
like) , a single source Si is associated with the 
surface sample dSi. As will be seen later, in the case 



L 
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of the estimation of a vector quantity in a point of 
space M, three sources will rather be assigned to each 
surface sample dSi. 

5 Preferably, the hemisphere HEMi is constructed as 
described hereinbelow. During the aforesaid meshing 
step a) , the surface area of the obstacle OBS is on the 
one hand evaluated, and, on the other hand, a number of 
surface samples dSi is chosen according to the desired 
10 accuracy of estimation of the physical quantity at the 
point M. Thus, the surface area of a sample dSi is 
given by Sq/N, where So corresponds to the total surface 
area of the obstacle and N corresponds to the chosen 
number of surface samples dSi. 

15 

The hemisphere HEMi has the same surface area as the 
sample dSi. Thus, the radius Ri of the hemisphere is 
deduced from the expression: 




20 

Each mesh cell represented by a surface sample dSi 
exhibits, in the example described, a parallelogram 
shape, whose centre Pi corresponds to the point of 

25 intersection of the diagonals of this parallelogram. 
The hemisphere HEMi is tangential to the surface sample 
dSi at this point Pi. Of course, the mesh cells may be 
of different shape, triangular or other. It is 
indicated generally that the point Pi corresponds to 

30 the barycentre of the mesh cell. 

The position of the source Si (situated at the centre 
of the hemisphere HEMi) is thus defined. The distance 
separating the source Si from the point of contact Pi 
35 corresponds to the radius Ri of the hemisphere HEMi and 
the straight line which passes through the points Pi 
and Si is orthogonal to the mesh cell dSi . 



J-. 
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In the example represented in Figure lA, the surface of 
a radiating element ER, corresponding for example to a 
wave generator, is furthermore meshed. With each mesh 
cell of the surface of the radiating element is 
5 associated a surface sample dS'i, as will be seen 
later. 



The matrix system that is formulated in the aforesaid 
step b) corresponds to: 

10 



^V(M|) ^ 






V(M2) 


= F X 




.V(Mk), 







[1] 



where : 

- the coefficients Vj (with 2^ --w N) of the 
15 first column matrix correspond to values each 

associated with a source, such as an electric charge 
(in the case of the estimation of an electric 
potential) , or to a magnetic flux (in the case of the 
estimation of a magnetic potential) , or else to a 
20 sound speed (in the case of the estimation of an 

acoustic pressure related to the propagation of a 
sound wave) ; 

- the coefficients V(Mi) (with i=l, 2, N) of the 
second column matrix each correspond to a value of 

25 the physical quantity (an electric or magnetic 

potential or a pressure) to be estimated at a point 
iyii of space) ; 

- the interaction matrix F comprises coefficients Ci,j, 
the general expression for which is given by: 

30 

Cu=f{MiSj) [2] 



It is thus understood that the coefficients of the 
matrix F are interaction coefficients which depend on 
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the distance separating each 
a source Sj associated with a 



point of the space 
mesh cell dS j . 



from 



10 



In the case of the propagation of an electric wave, the 
coefficients Ci,j, Vj and V{Mi), respectively of the 
interaction matrix of the first and of the second 
column matrix, are given by: 



c ^ 



2jtSQ MjSj 

Vj =qj [3] 
V(Mi) = Ui 



where : 

- So is a dielectric constant, 

- M^Sj is a distance measured as an algebraic value, 

- qj corresponds to an electric charge characterizing a 
15 source Sj, and 

- Ui corresponds to an electric potential at the point 
Mi. 

In the case of the propagation of a magnetic wave, the 
20 expression for these coefficients is as follows: 



iTTfiQ MjSj 

Vj =<Pj [4] 

V(Mi) = ^i 



where : 

- |io corresponds to the magnetic permeability of the 
25 medium where the point Mi is situated, 

- (Pj corresponds to the magnetic flux associated with 
the source Sj,- 
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- 9i corresponds to the magnetic potential at the point 
Mi. 



Within the framework of the propagation of an 
5 ultrasound wave, these coefficients are given by: 



c>.j = - 



(ik.MiSj) 



MiS 



dSi 



Vj = Vj 



[5] 



V(Mi) = p| 

in which: 

- ±2=-!. 

10 - (D is the angular frequency of the sound wave; 

- p is the density of the medium in which the point Mi 
is situated; 

- the vector Vj corresponds to the sound speed 
emanating from the source Sj; 

15 - k corresponds to the wave vector of the sound wave; 
and 

- Pi corresponds to the acoustic pressure generated by 
the propagation of the ultrasound wave at the point 
Mi. 

20 

In the expression for the coefficients ci,j, the term 
dSj corresponds to the surface area of the sample 
associated with the source Sj . Preferably, the meshing 
of a surface within the meaning of step a) of the 
25 method according to the invention is chosen in such a 
way that each mesh cell comprises one and the same 
surface area dS=dSi=dS2= . - . =dSj . 



It is noted in particular in the expression for the 
30 coefficients Ci,j that they depend on the scalar product 
of the wave vector and the vector MjSj . Thus, for 

ultrasound waves, account is taken of a phase shift 
between the paths which join each source Sj to a point 
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of three-dimensional space this phase shift being 

due to a difference in journey length between the rays 
leaving each source and arriving at the point M (as 
shown in Figure 4A) . In particular, it will be 
5 understood that the angle of incidence of such a ray is 
taken into account in the expression for the 
coefficients of the interaction matrix F. 



Of course, within the framework of the propagation of 
10 an electromagnetic wave of high frequency, hence of 
short wavelength, which differs from the electrostatic 
or magnetostatic framework hereinabove, account is 
taken of the propagation term exp(i^.f) in the 

expression for the interaction matrix, with respect to 
15 the geometry of the problem to be solved, as within the 
framework of the propagation of an ultrasound wave 
hereinabove (relation [5] ) . 



20 



25 



Thus, the matrix system of equation [1] makes it 
possible to estimate, on the basis of an interaction 
matrix F and of a vector comprising the values Vj 
associated with the sources Sj, the coefficients of a 
vector (column matrix) comprising the values of 
physical quantity V(Mi) at the points of space Mi. 



To determine the values of the 
following matrix system is applied: 



sources 



the 



rv(Pi) ^ 

V(P2) 



V(PN)j 



F X 



V2 



[6] 



30 



where : 

- the coefficients of the interaction matrix F are 
expressed by with Ci,j=f (PiSj) 
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- the indices i and j correspond respectively to the i^^ 
row and the j*^^ column of the interaction matrix F. 
This interaction matrix comprises, for the 
determination of the values associated with the 

5 sources Vj, N rows and N columns, recalling that N is 
the total number of mesh cells on the surface of the 
obstacle; 

- the points Pi correspond to the apex of the 
hemispheres HEMi of Figure IB. 

10 

The implementation of step c) of the method within the 
meaning of the present invention corresponds to 
computing a boundary condition for the points Pi, of 
known properties, as will be seen later. 

15 

The matrix system of equation [6] then becomes: 



V(P2) 

Iv(Pn)J 



t73 



where : 

20 - F""*" corresponds to the inverse of the interaction 
matrix F; and 

- the values V(Pi) are predetermined, as a function of 
the aforesaid boundary conditions. 

25 The source values Vj are thus determined. 



On the basis of the estimation of these source values 
Vj, it is possible to compute the scalar physical 
quantity at any point M of three-dimensional space, on 
30 the basis of the relation: 



[8] 
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To obtain this expression for the scalar quantity V(M), 
the interaction matrix F may comprise just one row of 
coefficients Cj, with: 
5 Cj = f(MSj), 

but always comprises N columns. 

Referring again to Figure lA, it will be understood 
that the surface of the obstacle OBS receiving the wave 
10 emitted by the radiating element ER acts, itself, as an 
active surface re-emitting a secondary wave (for 
example by reflection) . Each source Si represents a 
contribution to the emission of this secondary wave, 

15 Furthermore, to take account both of the presence of 
the main wave and of the presence of the secondary wave 
at the point M, the contribution of the main wave and 
the contribution of the secondary wave at the point M 
are estimated via the matrix system: 

20 









fv'i ^ 


V(M) « Fx 




+ Fx 


v'2 

• 











where : 

- F' is the matrix of interaction between the surface 
25 of the radiating element and the point M; 

- v*j (with j=l, 2, 3, N') is the value of the 
sources allocated to each surface sample dS ' j of the 
radiating element, N' being the chosen total number 
of mesh cells of the active surface of the radiating 

30 element ER, 

The coefficients of the matrix F' are again dependent 
on the distance MS'j, where S'j are the sources assigned 
to each surface sample dS'j of the radiating element. 
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According to an advantageous characteristic, the values 
of the sources of the obstacle Vj are determined as a 
function of the values of the sources of the radiating 
5 element v'j, and which are themselves computed as will 
be seen later with reference to Figures 4A, 4B, 5A and 
5B. 

Reference is now made to Figure 2A, in which three 
10 sources are assigned to each surface sample dSi, with a 
view to estimating a vector physical quantity V(M), at 
a point M of three-dimensional space. 

It will be understood that to estimate the vector 
15 quantity, via its three coordinates in space x, y and 
z, the number of equations to be solved with respect to 
the previous matrix system must be multiplied by three. 
Thus, the matrix F"""^ of relation [7] must comprise three 
times as many rows as before. The interaction matrix F 
20 must, itself, comprise three times as many columns as 
before and, accordingly, three sources are 
advantageously envisaged per mesh cell when dealing 
with the determination of the coordinates in three- 
dimensional space of a vector V(iyi). 

25 

Referring to Figure 2B, the three sources SAi, SBi, SCi, 
allocated to a surface sample dSi have respective 
positions determined as is indicated hereinbelow. Such 
as represented in Figure 2B, the three sources SAi, SBi, 

30 SCi are coplanar and the plane which includes these 
three sources also includes the base of the hemisphere 
HEMi. The hemisphere HEMi is constructed as indicated 
hereinabove (with the same surface area as the surface 
area of the mesh cell) , with however the centre of the 

35 hemisphere which corresponds here to the barycentre of 
the three sources SAi, SBi and SCi. 
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The "centre of the hemisphere" is understood to mean 
the centre of the disk which constitutes the base of 
the hemisphere. 

5 The three sources which are allocated to the surface 
sample dSi are placed at the vertices of an equilateral 
triangle whose barycentre Gi corresponds to the centre 
of the hemisphere. Preferably, each source SAi, SBi and 
SCi is placed in the middle of a radius Ri of the 
10 hemisphere. Thus, the straight lines which connect the 
barycentre Gi to each source are angularly separated by 
120° . 

Referring to Figure 2C, the angular orientation of the 
triangles formed by the triplets of sources is chosen 
randomly, from one surface sample to another* 
Advantageously, overperiodicity artefacts, which could 
result from the choosing of one and the same angular 
orientation of these triangles, are thus avoided in the 
estimation of the vector quantity at the point M. 

With reference to the various wave types indicated 
previously, the vector quantity V(M) to be estimated 
may be : 

25 - an electric field, within the framework of the 
propagation of an electric wave; 

- a magnetic field, within the framework of the 
propagation of a magnetic wave; and 

- a speed of sound at the point M, within the framework 
30 of the propagation of ultrasound waves. 



To determine the values 
SAi, SBi, SCi, the matrix 
to the following relation: 



associated with each source 
system is fomulated according 
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Vx(P2) 




vA2 


Vx(Pn) 




vAn 


Vy(Pl) 




vBi 


Vy(P2) 




vB2 


Vy(PN) 




vBn 


VzCPi) 




vCi 


Vz(P2) 




vC2 


^V^(Pn); 







[11] 



It is noted, in particular, that the interaction matrix 
FV is of dimensions 3Nx3N, where N is the total number 
of surface samples. The interaction matrix is expressed 
here through the relation: 



N 



N 



N 



Njci(iJ) C|(i,j) C£(i,j) 

N^D^cij) c^a,]) cy(ij) 

nHoj) C|(i,j) C^(i,j) 



[12] 



The coefficients of this matrix are expressed through: 



CS0J) = /«k/5.5^y)] [13] 

10 

with a = A, B, C 

i = 1, 2, . . . , N 
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j = 1, 2, . . . , N 
u = X, y, z. 



and are also dependent on a distance separating the 
5 point of contact Pi from one of the sources Saj (a = A, 
B or C) of a triplet associated with a surface sample 
dSj. 

By inverting the interaction matrix Fv t the values vaj 
10 associated with each source Saj are thus determined by 
applying boundary conditions on the values of • the 
vector V at the points Pi. These boundary conditions 
impose a value of the vector V, according to its three 
coordinates Vx(Pi), Vy(Pi) and Vz(Pi)- 

15 

Once these source values vaj have thus been determined, 
the expression for the vector V at any point M of space 
can easily be computed through the relation: 



20 



V(M) ^ Vx(M)x + Vy(M)y + V2(M)z 



[14] 



X, y and z correspond to unit vectors plotted along 
the axes x, y and z of three-dimensional space. 
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Thus, the interaction matrix Fv / when it is applied to 
any point M of space, ultimately comprises only three 
rows each associated with a coordinate of space x, y or 
z • 

5 

For various types of waves, the values of the sources 
vaj are, as before, an electric charge in respect of an 
electric wave, a magnetic flux in respect of a magnetic 
wave, a speed of sound in respect of an ultrasound 
10 wave. 

More precisely, the coefficients of the interaction 
matrix Fv are determined from the above relations [3] , 

[4] and [5], specifying however that: 

15 

V(M) = -gmd[V(M)] [15] 

V(M) being the scalar quantity computed previously 
through equation [8] . 

20 Thus, for the estimation of a vector quantity V at the 
point M and for the wave types mentioned above by way 
of example (electric, magnetic and ultrasound) , the 
coefficients of the interaction matrix Fv are inversely 

proportional to the square of a distance separating 
25 each source from the point M, while for the estimation 
of a scalar quantity V at a point M of space, the 
coefficients of the interaction matrix F are simply 
inversely proportional to this distance. Each distance 
involves one of the sources of a triplet of a surface 
30 sample and a point M of space. The interaction matrix 
Fv then comprises 3N columns when taking three sources 

per surface sample, while the interaction matrix F for 
the estimation of the scalar quantity comprised only N 
columns since just one source per surface sample was 
35 necessary. 
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More generally, one source per sample is allocated when 
boundary conditions are known for a scalar quantity and 
three sources per sample are allocated when boundary 
conditions are in fact known for a vector quantity. 

5 

Reference is now made to Figure 3A to describe, by way 
of illustration, an application of the method according 
to the invention to the estimation of an electric 
potential at a point M of three-dimensional space, 

10 situated between two plates of a capacitor. The plates 
of this capacitor are brought to respective potentials 
VI and V2 . The implementation of step a) consists 
firstly in meshing the respective surfaces of the two 
plates. In the example represented in Figure 3A, only 

15 two mesh cells have been represented for each plate, 
simply by way of illustration. 

The application of step b) consists in formulating the 
matrix system involving the interaction matrix F and 
20 the column vector comprising the values of the sources 
Si to S4. Multiplication of these two matrices makes it 
possible to obtain a column vector comprising the 
values of the potential at one or more points M of 
space . 

25 

The implementation of step c) of the method according 
to the invention consists in applying the matrix system 
to the points of contact of the hemispheres Pi to P4, of 
each surface sample dSi to dS4. This results in the 
30 following relation: 



'V(P,)) 

V(P2) 
V(P3) 

.V(P4)J 



2neQ 



PlSi 
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1 
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P1S3 

1 

P2S3 

1 

P3S3 

1 



1 

P1S4 

1 

P2S4 

1 

P3S4 

1 



^^P4Si P4S2 P4S3 .P4S4^ 



V2 
V3 
^^4, 



[16] 
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with 



V(Pi) = V{P2) = Vi 

V(P3) = V(P4) = V2 

Vi = qi, V2 = q2/ V3 = q3, V4 = q4 



5 



Here, the boundary condition prescribes that the value 
of the potential at the points of contact Pi and P2 
should correspond to the potential VI of the first 
plate. Likewise, the electric potential at the points 

10 of "contact P3 and P4 should correspond • to the electric 
potential of the second plate V2 . By inverting the 
interaction matrix applied to the point of contact Pi, 
the values of the sources vi which correspond, as 
expressed in relation [16], to electric charges qi are 

15 determined. 

The coefficients of the interaction matrix are 

27l8oPiSj 

known perfectly, since the positions of the sources Sj 
and the positions of the points of contact Pi are 
20 determined beforehand, as is represented in Figure IB. 

The expression for the electric potential V(M) at the 
point M between the two plates is ultimately given by 
the expression: 



Reference is now made to Figure 3B in which the same 
plates have been represented as in Figure 3A, with 
30 substantially the same mesh, but with the aim, here, of 
estimating a vector quantity corresponding to the 
electric field E(M) , at the point M of three- 
dimensional space . 



25 



V(M) = 




C17] 
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Relations [11] to [15] can be applied to estimate the 
value of the electric field at the point M, with, in 
relation [13] : 



[18] 



with 
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= A, 


B, 


C 
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= X, 
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However, the values of the electric field at the point 
of contact Pi remain to be determined in relation [11]* 

15 

A predetermined general law for the behaviour of the 
field (in reflection, in transmission or other) at the 
level of the surface of the obstacle (plates in the 
example of the aforesaid capacitor) is then introduced 
20 to ascertain the values of the sources vaj . 

For example, if the electric wave is reflected totally 
by the surface of an obstacle (for example one of the 
two plates) , the electric field at a point of contact 

25 Pi is normal to the surface dSi and its components Ex 
and Ey are zero* By way of illustration, if the surface 
of the plate was represented only by a single surface 
sample with three sources, the values of its sources 
vA, vB and vC would all be mutually equal to one and 

30 the same value +q. 

On the other hand, if the coefficient of reflection is 
practically zero at the surface dSi, the component of 
the electric field Ez at the point Pi is zero, this 
35 corresponding indeed to the case where the field is 
substantially tangential to the surface dSi. Thus> by 
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way of illustration, if the surface of the plate was 
represented only by a single surface sample with three 
sources, the values of its sources vA, vB and vC would 
be, for example, +q, +q and -2q respectively. For 
5 example, within the framework of the propagation of a 
magnetic wave, if the surface of an eddy current sensor 
(with a zero normal component of the magnetic field) 
was represented by a single surface sample, the 
magnetic fluxes of the three sources associated with 
10 this surface sample would be +(p, +cp and -2q>. 

It is thus understood that with the three sources per 
sample dSi, it is possible to define, for example as a 
function of the weighting of each source, any 
15 orientation of the field at the surface of the 
obstacle . 

Of course, this approach assumes that the coefficient 
of reflection R of an obstacle is known beforehand. In 

20 particular, it may be advantageous to compare a 
simulation and an experimental measurement so as to 
detect, at the surface of an obstacle, inhomogeneities 
or impurities which correspond to points of the surface 
of this obstacle which do not satisfy the values of the 

25 coefficient of reflection R that are prescribed at each 
predetermined point Pi of the obstacle. 

A predetermined value of the reflection coefficient can 
thus be assigned to each point Pi of the surface of the 

30 obstacle ♦ Accordingly, a matrix R which is 
representative of the reflection coefficient at each 
point Pi is introduced. For an interaction between a 
radiating element and an obstacle, it is thus possible 
to express the matrix .system of relation [9] in another 

35 way, that is to say by giving a single expression for 
all the sources of the system (both of the obstacle and 
of the radiating element) , as indicated hereinbelow. 

In what follows, it is indicated that: 
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- F(P) is the interaction matrix of the obstacle OBS, 
applied to the points Pi of the surface of the 
obstacle OBS; 

- F(P') is the interaction matrix of the obstacle OBS, 
5 applied to the points P'l of the surface of the 

radiating element ER; 

- F' (P) corresponds to the interaction matrix of the 
radiating element ER, applied to the points Pi of the 
surface of the obstacle OBS; 

10 - F'(P') corresponds to the interaction matrix of the 
radiating element ER, applied to the points Pi of the 
surface of the radiating element ER; 

- V ' corresponds to the column vector comprising the 
values of the sources S'i of the radiating element 

15 ER; and 

- V corresponds to the column vector comprising the 
values of the sources Si of the obstacle OBS. 

On an obstacle, the contribution of the wave emitted by 
20 the radiating element ER is expressed by: 



V'(P) = F(P).v' [19] 

The contribution of the secondary wave returned by the 

obstacle OBS is expressed, by definition, by the 
25 relation: 



V(P) = F(P).v C20] 

Now, in the example represented in Figure 4A, the 
secondary wave corresponds simply to a reflection of 
30 the main wave. This is expressed by the relation: 



V(P) = R V'(P) 



[213 
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iO 



15 



20 



25 



where R corresponds to a reflection matrix each 
coefficient of which represents the contribution to the 
emission, by reflection, of the secondary wave, by each 
source Si (or Sa±, within the framework of estimation of 
a vector quantity) of the obstacle OBS. 

From the three relations [19] , [20] and [21] is deduced 
the expression for the column vector v comprising the 
values of the sources of the obstacle, on the basis of 
the column vector v comprising the values of the 
sources of the radiating element, through the relation: 



Additionally, for fine estimation of the scalar or 
vector quantities at the point M, in particular to take 
account of multiple reflections, it is advantageous to 
take account of the contribution of the radiation by 
the obstacle, at the level of the surface of the 
radiating element ER. Accordingly, account is taken, in 
the estimation of the boundary conditions at the 
surface of the radiating element ER (at the points P'l) 
of the contribution of the radiation of the sources S'l 
of the radiating element and of the contribution of the 
emission of the secondary wave by the sources Si of the 
obstacle, through the relation: 



The source values S'i of the radiating element ER can 

thus be tailored, by virtue of relation [23] , by taking 

account of the reflection of the obstacle OBS, 
according to the following relation: 




Vt(F) = F(F) v-f F(F) 



[23] 
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Vt(P*) - ^(P*) AF(P) .R ,[r(P) ] + F'(P') ],v' [10] 

Thus, boundary conditions are simply prescribed for the 
radiating element, so as to deduce therefrom the values 
of the sources v'l. In practice, one will preferably 
5 proceed as follows: 

- after meshing the surfaces, one determines the 
position of the points Pi and P'l and of the sources Si 
and S'i; 

- as a function of the type of wave involved, one 
10 determines the coefficients of the matrices F(P), 

F' (P) , F(P' ) and F' (P' ) ; 

- as a function of a law of reflection of the obstacle, 
one determines the coefficients of the reflection 
matrix as in the example given later in respect of an 

15 ultrasound wave; 

- as a function of boundary conditions on the radiating 
element (the behaviour of which is generally known for 
a given problem) , one determines the values of the 
vector Vt(PM at the points P'i of the surface of the 

20 radiating element and one deduces therefrom the values 
of the sources S'l of the radiating element by 
inverting relation [10] ; 

- one also deduces therefrom the values of the sources 
Si of the obstacle by applying relation [22]; 

25 - once the values of all the sources Si and S'i have 
been determined, one can apply the matrix system given 
by relation [9] to any point M of space, by applying 
the interaction matrices F and F' (involving the 
position of the point M and the positions of the 

30 respective sources Si and S'i) to this point M. 

Referring again to Figure 4A, it is considered that the 
obstacle OBS represents simply an interface between two 
media Ml and iyi2, thus forming a dioptric member which 
35 may be plane, such as represented in the example in 
Figure 4A, but also curved or of any general shape. The 
reflection coefficients Ri associated with each point Pi 



- 29 - 

depend, within the framework of the propagation of an 
ultrasound or electromagnetic wave of high frequency, 
on the angle of incidence pi of the ray emanating from 
the source Si, at the point of three-dimensional space 
5 M. 



For an ultrasound wave, the expression for the 
coefficients of reflection Ri is given by: 



p2<it COSfii - piC| 




1/2 


P2C2 COS A + piCi 


. cf cf 


1/2 



[24] 



10 



where: 

- ci is the speed of sound in medium Mi; 

- C2 is the speed of sound in medium M2; 
15 - pi is the density of medium Mi; 

- P2 is the density of medium M2. 

In this expression [24], the term cosPi may simply be 
estimated as a function of the coordinates in space of 
the point M and of the point representing the source 



Referring now to Figure 4B, the same estimation may be 
undertaken for a point M situated in the medium M2 . In 
this case, the wave that the point M receives is a wave 

25 transmitted by the obstacle OBS • In particular, it is 
noted that the sources of the radiating element ER are 
no longer active, owing to the occultation of the 
radiating element ER by the obstacle OBS. In 
transmission, the reasoning applies as above with a 

30 boundary condition prescribed at the points Pi through 
the values of the transmission coefficients Ti 
associated with each point Pi. Within the framework of 
the propagation of an ultrasound wave, each 
transmission coefficient Ti is given by the relation: 
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Ti = 



[25] 



1/2 



P2C2 cos/?i -I- pici cos^ 

L cf cf 



As indicated above^ the terms cos(3i may be determined as 
5 a function of the respective coordinates of the sources 
Si and of the point M. 

To estimate the values of sources Si of the obstacle 
OBS, relation [22] is applied while replacing, however, 
10 the reflection matrix R by the transmission matrix T: 



Within the framework of the propagation of an 
15 ultrasound wave, the coefficients of the matrices R and 
T are estimated for each source Si and for each point 
Pi. In particular, each coefficient Ti,j or Ri,j of the 
matrix T or of the matrix R (where i corresponds to the- 
i^^ row and j corresponds to the j^^ column) is 
20 expressed as a function of an angle Pij between a normal 
to the surface of the obstacle at the point Pi and a 
straight line passing through the point Pi and through 
a source Sj • It is thus possible to write, in a general 
manner, the two relations expressing the values of the 
25 coefficients of the matrices T and R by the following 
respective relations : 




[26] 



Tjj = ft (cosAj) 



[27] 



= ff (cos/ffy) 
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[283 



where ft is given by relation [25] and fr is given by 
relation [24] • 

5 

More generally, with reference to Figures 4A and 4B, it 
is indicated that, if the obstacle is considered to be 
a solid material representing a medium M2 distinct from 
a medium Ml in which the main wave propagated 
10 initially: 

- for a reflection of the main wave off the* obstacle in 
the guise of medium M2 (the surface of the obstacle 
forming a dioptric member between the media Ml and M2), 
the hemispheres HEMi are oriented outwards from the 

15 obstacle (Figure 4A) ; 

- for a transmission of the main wave within the 
obstacle, the hemispheres HEMi are oriented inwards 
into the obstacle (Figure 4B) . 

20 Reference is now made to Figure 5A to describe the case 
of a plane obstacle OBS of finite dimensions, excited 
by a radiating element ER, inclined by a predetermined 
angle with respect to the obstacle OBS. As indicated 
above, for an ultrasound wave, the inclination of the 

25 radiating element will be talcen into account to compute 
the contribution of the wave emitted by the radiating 
element at the point M. Additionally, in a particularly 
advantageous manner, a surface which encompasses the 
surface of the obstacle is meshed (Figure 5A) . For a 

30 slice of space that is delimited by the radiating 
element, on the one hand, and the obstacle, on the 
other hand (Figure 5A) , three types of sources can be 
considered : 

- the sources S'i of the radiating element ER, 

35 - sources SOi, which return the secondary wave, by 
reflection from the obstacle OBS, as a function of a 
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certain coefficient of reflection R of the obstacle; 
and 

- sources SSi, which do not return any secondary wave 
and to which a zero coefficient of reflection may be 

5 allocated if the obstacle separates two media of 
identical indices. In this case, these sources SSi 
are considered to be "off" in the aforesaid slice of 
space and are not taken into account in the 
computations of the physical quantity at the point M 
10 of Figure 5A, On the other hand, these sources SSi 

may be active through reflection of the main wave if 
the obstacle OBS separates two media of different 
indices, 

15 Additionally, to estimate the scalar or vector 
quantities associated with a point M of a half-space 
delimited by the surface encompassing the obstacle OBS 
(on the right of Figure 5B) , one considers: 

- the sources SO'i of the obstacle, which are active 
20 through transmission of the main wave, and 

- the sources SSi, to which one now assigns a 
coefficient of transmission equal to 1 if the 
obstacle separates two media of like indices. These 
sources SSi behave ultimately (to within angles of 

25 incidence) like the sources S'i of the radiating 

element ER. 

The sources S'i of the radiating element may then be 
"off" for the computation of the physical quantities in 
this half-space, 

30 

To compute the values v' of the sources S'l of the 
radiating element (from which are deduced the values v 
of the sources of the obstacle according to relations 
[22] and [26]), boundary conditions will simply be 
35 applied to the points of the active surface of the 
radiating element ER. For example, for ultrasound wave 
propagation, it may be indicated that the sound 
velocities at the points of the surface of the 
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radiating element ER are perpendicular to this surface 
and have mutually equal moduli vq. 



Generally, it is indicated that the three-dimensional 
5 space may thus be divided up by interfaces delimiting 
media of distinct properties, each interface 
representing an obstacle within the meaning of the 
present invention. The physical quantities are then 
computed in each slice of space. For example, within 

10 the framework of the study of a heterostructure (for 
several interfaces), the above method may be applied in 
respect of successive slices of space by considering 
two interfaces: one representing a "radiating element", 
within the meaning of Figures 4A and 5A, for example 

15 through transmission of a received wave, and the other 
representing an obstacle receiving the transmitted 
wave. Advantageously, account is taken, for each slice 
of space^ of the contributions of all the interfaces, 
as expressed by relations [10] and [22] . 

20 

However, in a preferred practical embodiment, in 
particular in order to programme the simulation of an 
interaction, one will advantageously consider all the 
obstacles of the entire space around a point M and a 
25 condition will be prescribed regarding the position of 
the point M with respect to each source present in the 
space . 

Preferably, with reference to Figure 7A, the scalar 
30 product SM. f is tested at each iteration with respect 
to a source S, for example in the form: 



35 where r is the vector connecting the source S to the 
point of contact P of the half-sphere with the surface 
element dS considered, in the case where just one 
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source per hemisphere is envisaged. In the case where a 
triplet of sources SI, S2, S3 is in fact envisaged, the 
base of the vector f is preferably situated at the 
barycentre of the three sources SI, S2, S3. Moreover, 
5 in the case of three sources per surface sample, the 
computation of the scalar product concerns each source 
Si of the triplet SI, S2, S3. 

Typically, if cos9 is positive, one takes into account 
10 the contribution of this source S in the estimation of 
the interaction. 

On the other hand, if cosG is negative, a zero (scalar 
or vector) value is assigned to this source S in the 
15 estimation of the interaction. 

As a variant of the computation of the scalar product 
above, an "altitude" of the point M can be computed. In 
this case, the test pertains to a quantity of the type: 

20 




Of course, other types of test are possible. For 
example, in the case of a computation with respect to 
25 the angle 0, this angle can be chosen within a cone of 
chosen aperture, or otherwise. 

Ultimately, this approach advantageously makes it 
possible to systematize any configuration of the 
30 sources with respect to the observation point M, by 
simply introducing an extra step for testing, at each 
iteration on a source S, the position of this source S 
with respect to the point M, as indicated above. 

35 This approach proves to be particularly advantageous 
for surfaces to be meshed which are relatively complex, 
in particular when the observation point M is liable to 
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be situated in a shadow zone with respect to certain 
sources, as shown in Figure 7B. In this Figure 7B, the 
half-spheres associated with the sources in the shadow 
zone of the observation point M and for which, 
5 consequently, the contribution is fixed as being zero 
in the estimated interaction have been represented by 
dashes. Of course, in the case where a source is 
screened by a sampled surface, although the scalar 
product SM. f associated with this source remains 
10 po-sitive, a second test determines whether - the vector 
SM does or not does cross a sampled surface. If it 
does, this source is considered to be inactive 
specifically for the region of the point M. 

15 Thus, in more general terms, the method within the 
meaning of the invention preferably envisages at least 
one extra step, for each surface sample, for testing 
the value of a scalar product of: 

20 - a first vector f normal to the surface sample and 
directed towards the apex P of the hemisphere, such as 
represented in Figure 7A, and 

- a second vector SM drawn between a source S 
25 associated with this hemisphere and the point M which 
is situated in the region of observation, 
while distinguishing, in particular: 

- the case where this scalar product is less than a 
predetermined threshold and the contribution of 

30 this source is not taken into account, and 

- the case where this scalar product is greater than 
a predetermined threshold and the contribution of 
this source is actually taken into account. 

35 In the example above where the angle 0 between these 
two vectors is considered, the aforesaid predetermined 
threshold is of course the value zero and one simply 
distinguishes between the cases where the scalar 
product is positive or negative. 
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Of course, this choice is not limiting so that, for a 
heterostructure with several parallel dioptric members, 
it will again be possible to consider, advantageously, 
5 successive half-spaces, as described above with 
reference to Figures 4A and 4B. 

The simulation of Figure 5C corresponds, for an 
ultrasound wave, to the situation of Figures 5A and 5B 
10 by taking accouTit: 

- of the contribution of the emission of the main wave 
by the radiating element ER; 

- of the contribution of the reflection of this main 
wave by the obstacle; and 

15 - of the contribution of the transmission of the main 
wave by the obstacle. 

The level lines of Figure 5C correspond to various 
magnitudes of acoustic pressure. The radiating element 

20 ER is placed 10 mm from the obstacle OBS and inclined 
by 20*^ with respect to the latter. Interference fringes 
are noted in particular in a zone between the obstacle 
OBS and the radiating element ER. Such a simulation may 
advantageously indicate an ideal position of an 

25 ultrasound sensor. These ultrasound sensors customarily 
comprise a transducer as active radiating element and a 
detector for measuring the ultrasound waves received. 
The simulation of Figure 5C may thus moreover indicate 
the ideal shape of an ultrasound sensor, depending on 

30 the desired applications, for a given shape of 
obstacle . 

The simulation of Figure 5C has been performed by 
virtue of a matrix computation programmed with the aid 
35 of the MATLAB© computation software. The chosen total 
number of mesh cells for the obstacle and for the 
radiating element (here, a few hundred in all) is then 
optimized: 
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- on the one hand, to limit the duration of the 
computations; and 

- on the other hand, so that the size of the mesh cells 
remains less than half a wavelength, so as to satisfy 

5 the Rayleigh criterion. 

It is indicated however that, as the elements to be 
meshed in the implementation of the method according to 
the invention are simply surfaces, the computation 
10 • tdiues are not nearly as long as those required in the 
implementation of a method of computation of "finite 
element" type. 

The present invention can thus be realized by the 
15 implementation of a succession of instructions of a 
computer program product stored in the memory of a hard 
disk or on a removable support and running as follows: 

- choice of a meshing stepsize in particular as a 
function of the wavelength of the main wave; 

20 - determination of the coordinates of the sources Si 
and/or S'l and of the points of contact Pi and/or P'i; 

- choice of a type of wave involved and computation of 
the coefficients of the interaction matrices applied 
to the points Pi and/or P'i through the implementation 

25 of matrix computation software; 

- choice of a law of reflection and/or of transmission 
of the surface of the obstacle and computation of the 
coefficients of the reflection and/or transmission 
matrices ; 

30 - computation of the values of the sources S± and/or 
S'i; and 

- computation of the scalar or vector quantities at 
every point of three-dimensional space. 



35 In this regard, the present invention is also aimed at 
such a computer program product, stored in a central 
unit memory or on a removable support able to cooperate 
with a reader of this central unit, and comprising in 
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particular instructions for implementing the method 
according to the invention. 



Of course, the present invention is not limited to the 
5 embodiment described hereinabove by way of example; it 
extends to other variants. 

Thus, it will be understood that, even if, in the 
figures discussed above, both the surface of an 

10 obstacle and the surface of a radiating element are 
represented, the present invention applies also to the 
estimation of physical quantities within the framework 
of a wave interacting with an obstacle and emitted in 
the far field. In this context, it is not necessary to 

15 demarcate the surface of a radiating element to be 
meshed and relations [8] and [14] above suffice to 
determine the interaction between this wave and the 
obstacle . 

20 Equations making it possible to compute the scalar or 
vector quantities at a point M of space, for 
electromagnetic waves, or acoustic waves, have been 
indicated above. Of course, these quantities may be 
estimated for other types of waves, in particular for 

25 thermal waves, electromagnetic waves involving 
radiof requency antennas, or others. 

Of course, the present invention is not limited to an 
application to nondestructive testing, but to any type 
30 of application, in particular in medical imaging, for 
example for the study of microsystems employing 
acoustic microscopy with movable mirrors. 

Interactions between a wave and a single obstacle have 
35 been described above. Of course, the present invention 
applies to an interaction with several obstacles. 
Accordingly, it is simply necessary to mesh the 
surfaces of these obstacles and to add up their 
contribution for the estimation of a vector or scalar 



quantity at any point of space. Likewise, as indicated 
above, the surface of the obstacle OBS may be plane, or 
else curved, or else of any complex shape. 

Thus, within the framework of a wave interacting with 
several obstacles in space, a simulation equivalent to 
that represented in Figure 5C would make it possible to 
position sensors and/or radiating elements as a 
function of the configuration of these obstacles, in 
pa-rticular for an application to the 'determination of 
the position of loudspeakers in a partitioned cabin, 
such as a motor vehicle cabin. 

The three-dimensional space may be divided up into a 
plurality of regions, as described above with reference 
to Figures 4A, 4B, 5A and 53. However, in order for a 
surface of one of the said regions to be considered to 
be an obstacle of a main wave, the said obstacle 
becoming active by emission of a secondary wave, the 
angle of incidence of the main wave on this surface 
must preferably remain less than or equal to 90°. 



